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Observations of galaxy clusters show that the intracluster medium (ICM) is likely to be turbulent 
and is certainly magnetized. The properties of this magnetized turbulence are determined both by 
fundamental nonlinear magnetohydrodynamic interactions and by the plasma physics of the ICM, 
which has very low collisionality. Cluster plasma threaded by weak magnetic fields is subject to 
firehose and mirror instabilities. These saturate and produce fluctuations at the ion gyroscale, which 
can scatter particles, increasing the effective collision rate and, therefore, the effective Reynolds 
number of the ICM. A simple way to model this effect is proposed. The model yields a self- 
accelerating fluctuation dynamo whereby the field grows explosively fast, reaching the observed, 
dynamically important, field strength in a fraction of the cluster lifetime independent of the exact 
strength of the seed field. It is suggested that the saturated state of the cluster turbulence is 
a combination of the conventional isotropic magnetohydrodynamic turbulence, characterized by 
folded, direction-reversing magnetic fields and an Alfven-wave cascade at collisionless scales. An 
argument is proposed to constrain the reversal scale of the folded field. The picture that emerges 
appears to be in qualitative agreement with observations of magnetic fields in clusters. 



I. INTRODUCTION 

Clusters of galaxies are vast and varied objects that 
have long attracted the attention of both observers (in 
recent decades spectacularly aided by X-ray and radio 
telescopes) and theoreticians. The observed properties 
of clusters have proved far from easy to explain as new 
data has confounded many old theories. 

The overall budget of the cluster constituents is 
roughly as follows: ~ 75% of cluster mass is dark mat- 
ter, whose sole function is assumed to be to provide the 
gravitational well, ~ 20% of cluster mass is the diffuse X- 
ray-emitting plasma (the intracluster medium, or ICM), 
while the galaxies have an all but negligible mass. The 
plasma that makes up the ICM is made of hot and ten- 
uous ionized hydrogen: temperatures are in the range 
of 1 — 10 keV, number densities ~ 10 _1 — 10~ 3 cm -3 , 
with colder, denser material found in the cool cores and 
hotter, more diffuse one in the outlying regions^ 

Necessarily, the first observations and the first the- 
oretical models of clusters concerned what one might 
call large-scale features, such as the overall profiles of 
mass and temperature, the structure formation and the 
role of the central objects. As observations increased 
in accuracy and resolution, the ICM was revealed to 
be much richer than simply a dull cloud of X-ray glow 
smoothly petering out with distance from the center. 
A great panoply of features has been detected: bub- 
bles, filaments, ripples, edges, shocks, sound waves, etc., 
as well as very chaotic density, temperature and abun- 
dance distributions It is particularly the presence 
of chaotic fields and the evidence that this chaos exists 
in a range of scales^! that makes one expect that ICM, 
like so many other astrophysical plasmas, is in a turbulent 
state. It is essential to know the properties of this turbu- 
lence in order to predict the current and future statistical 



measurements of the plasma and magnetic fields (spectra, 
correlation and distribution functions, etc.) and to model 
correctly the transport processes in the ICM that deter- 
mine, for example, the overall temperature profilesAS 

We shall assume that the physics of small scales is, 
at least to some degree, independent of large-scale cir- 
cumstances and that we can, therefore, gain some use- 
ful understanding of the turbulence in clusters by ignor- 
ing large-scale features and considering a homogeneous 
subvolume of the ICM. In what follows, after reviewing 
briefly what is known about fluid motions (Sec. [HJ and 
magnetic fields (Sec. llUfl in clusters, we describe, mostly 
in qualitative terms, what we consider to be the essential 
aspects of the small-scale physics of the ICM (for plasma 
physics, see Sec. I1V|) . This will lead us to a tentative 
overall picture of the structure of the ICM turbulence 
(Sec. IVI|) as well as of the origin of its magnetic compo- 
nent (Sec.EJ. We must emphasize that the current state 
of the debates on the nature of turbulence in clusters 
(and, indeed, on its very existence^) is such that even 
a fundamental, conceptual view of the problem has not 
yet been agreed, and, therefore, a quantitative theory — 
even an incomplete one such as exists for hydrodynamic 
turbulence — remains a matter of future work. 



II. TURBULENCE 

The failure of one of the instruments on the ASTRO- 
E2 satellite has set off the planned direct detection of 
cluster turbulenceiii^ into the (probably not very dis- 
tant) future. However, indirect evidence of turbulent gas 
motions does exist: three recent examples are the broad 
spectrum of pressure fluctuations measured in the Coma 
cluster—, detection of subsonic gas motions in the core of 
the Perseus cluster^ and, again for the Perseus cluster, 
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TABLE I: Cluster Parameters 



Parameter 


Expression 


Cool cores" 


Hot ICM 


T 


observed 


3 x 10 7 K 


10 8 K 


n 


observed 


6 x 10~ 2 cm" 3 


10" 3 cm" 3 


«th,i 


(2T/m0 1/2 


700 km/s 


1300 km/s 


Vii 


1.5 nT- 3 ' 21 


5 x 10~ 13 s" 1 


2 x 10" 15 s~ 


Amfp 




0.05 kpc 


30 kpc 


/''II 


fth,iA m fp 


10 28 cm 2 /s 


10 31 cm 2 /s 


n 


3 x lo 13 ! 1-3 / 26 


200 cm 2 /s 


30 cm 2 /s 


u 


inferred 


250 km/s 


300 km/s 


L 


inferred 


10 kpc 


200 kpc 


L/U 


inferred 


4 x 10 7 yr 


7 x 10 8 yr 


Re 


uL/mi 


70 


2 


Rm 


ul/t, 


4 x 10 27 


6 x 10 29 


^visc 


(L/U)Re- 1/2 


5 x 10 6 yr 


5 x 10 8 yr 


^visc 


LRe" 3 / 4 


0.4 kpc 


100 kpc 


lies 


LRm- 1 / 2 


5000 km 


8000 km 




eB eq /crrii 


0.3 s" 1 


0.04 s" 1 




Vth 7 i /^z,cq 


3000 km 


30, 000 km 


B 


B e qPi,eq/A m fp 


5 x 10~ 17 G 


2 x 10" 19 G 


Bi 


Eq. COP 


3 x 10~ 14 G 


2 x 10" 17 G 


B 2 


Eq. nn e 


8 x 10~ 7 G 


2 x 10" 7 G 


-Bvisc 


B eq Re" 1/4 


9 x 10~ 6 G 


4 x 10" 6 G 


B COl 


(87rmin(7 2 /2) 1/2 


3 x 10~ 5 G 


4 x 10~ 6 G 






8 


20 


i± 


(B 2 /B cq )L 


0.2 kpc 


7 kpc 


Id 


observed 


1 kpc 


10 kpc 



"These numbers a re b ased on the parameters for the Hydra A 
cluster given in Ref. L'lt 

'in these expressions, n is in cm -3 , T in Kelvin. 

c We used a = 3/2 to get specific numbers, but the outcome is not 
very sensitive to the value of a. 



the broadening (assumed to be caused by turbulent dif- 
fusion) of abundance peaks associated with the brightest 
cluster galaxies^ These and other studies and models 
based on observational data appear to converge in ex- 
pecting turbulent flows with rms velocities in the range 
U ~ 10 2 — 10 3 km/s at the outer scales L ~ 10 2 kpc. 
The energy sources for this turbulence are probably the 
cluster and subcluster merger events and/or, especially 
for the turbulence in cool cores, the active galactic nu- 
clei (AGN). The aforementioned observational estimates 
of the strength and scale of the turbulence are in order- 
of-magnitude agreement with the outcomes of numerical 
simulations of cluster formationii*i£*i£i and of the buoy- 
ant rise of radio bubbles generated by the AGNpiLi^ Fur- 
ther discussion and references on the stirring mechanisms 
for cluster turbulence can be found in Rcfs. IH 13, EH 

While estimates of the turbulence parameters appear 
robust roughly to within an order of magnitude, a more 
quantitative set of numbers is elusive, partly because 
of the indirect and difficult nature of the observations, 



partly because the conditions vary both in different clus- 
ters and within each individual cluster. Here we shall 
adopt two fiducial sets of parameters: one for cool cores 
and one for the bulk of hot cluster plasma. These are 
given in Table [I] (along with some theoretical quantities 
that will arise in Sec. IVI and Sec. IV1[) . They will allow 
us to make estimates that will have the virtue of being 
consistent and systematic but must not be interpreted as 
precise quantitative predictions. They are representative 
of the range of conditions that can be present in clus- 
ters. The turbulence is assumed to be stirred at the outer 
scale L (with rms velocity U at this scale) and to have 
a Kolmogrov-type cascade below this scale. The small- 
scale cutoff is determined by the microphysical proper- 
ties of the cluster plasma. In Table HJ we give the value 
of the particle mean free path A m f p , which can be used 
in a naive estimate of viscosity: /iicM ~ Wth,iA m f p , where 
vth.i — (ZT/rrii) 1 / 2 is the ion thermal speed. We see 
that this gives fairly low values for the Reynolds number, 
Re ~ UL/hicm- R is this feature of the ICM that con- 
tinues to fuel doubts about its ability to support turbu- 
lence, at least in the strict, hydrodynamic high-Reynolds- 
number sensed However, one should be cognizant of the 
fact that whatever type of turbulence might exist in the 
cluster plasma, it is certainly not hydrodynamic, because 
this plasma is highly electrically conducting^ and mag- 
netized. The presence of the magnetic fields not only has 
a dynamical effect on the turbulence (due to the action 
of the Lorentz force, the medium acquires a certain elas- 
tic quality), but also changes the transport properties of 
the plasma itself: the viscosity, in particular, becomes 
strongly anisotropic^ These issues will constitute the 
main subject of this paper, but first let us briefly de- 
scribe what is known about magnetic fields in clusters. 



III. MAGNETIC FIELDS 

The first observed signature of cluster magnetic fields 
was the diffuse synchrotron radio emission in the Coma 
cluster detected in 1970i^ Starting from early 1990's, 
increasingly detailed measurements of the Faraday Ro- 
tation in the emission from intracluster radio sources 
have made possible quantitative estimates of the mag- 
netic field strength and scales in a large number of 
clusters Randomly tangled magnetic fields with 
rms strength of order B lms ~ 1 — 10 /iG are consistently 
found, with the fields in the cool cores of the cooling- 
flow clusters somewhat stronger than elsewhere. This is 
fairly close to the value _B cq that corresponds to magnetic 
energy equal to the energy of turbulent motions (see Ta- 
ble©. Thus, the magnetic field must be dynamically 
important. The estimates for the tangling scale Is of the 
field are usually arrived at by assuming that direction 
reversals along the line of sight (probed by the Faraday 
Rotation measure) can be described as a random walk 
with a single step size equal to Ib (the estimate of B rms 
is obtained in conjunction with this model). This gives 
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l B ~ 1- 10 kpc. 

The single-scale model is almost certainly not a correct 
description on any but a very rough level. Fortunately, 
much more detailed information on the spatial structure 
of the cluster fields is accessible. First, using certain sta- 
tistical assumptions (most importantly, isotropy), it is 
possible to compute magnetic-energy spectra from the 
maps of the Faraday Rotation measure associated with 
extended radio sources (the radio lobes of the jets emerg- 
ing from the AGN — these can be as large as ~ 10 2 kpc 
across) £21 This has been done most thoroughly for a ra- 
dio lobe located in the cool core of the Hydra A cluster^ 
The spectrum has a peak at k ~ 2 kpc -1 followed by 
what appears to be a power tail consistent with fc -5 / 3 
down to the resolution limit of k ~ 10 kpc -1 . The rms 
magnetic field strength is B ims = 7 ± 2 pG. 

The second source of information on the cluster field 
structure is the polarized synchrotron emission, which 
probes the magnetic field in the plane perpendicular to 
the line of sights Such data, while widely used for Galac- 
tic magnetic field studies^ has until recently not been 
available for clusters. This is now changing: the first 
analysis of polarized emission from a radio relic in the 
cluster A2256 reveals the presence of magnetic filaments 
with field reversals probably on ~ 20 kpc scale, which, 
however, is dangerously close to the resolution scaled 
This data is representative of the situation in the bulk 
of the ICM, rather than in the cores. Statistical analysis 
of such data will make possible quantitative diagnosis of 
the field structure and its dynamical rolei^I 

Thus, our knowledge of the magnetic fields in clusters, 
while far from perfect, is more direct and more detailed 
than that of the turbulent motions of the ICM. It is also 
due to improve dramatically with the arrival of new radio 
telescopes such as LOFAR and SKAil 3 - 



IV. PLASMA PHYSICS 

The key property of the ICM as plasma is that it is 
only weakly collisional and magnetized: given the ob- 
served values of the magnetic field, the ion gyroradius 
is pi ~ 10 4 km, which is much smaller than the mean 
free path. As p L <C A m f p already for dynamically very 
weak fields (B 3> Bq, see Table QJ), this is true both in 
the observed present state of the ICM and during most 
of its hypothetical past, when the magnetic field was be- 
ing amplified from some weak seed value. In a plasma 
with pi <§C A m f p , the equations for the flow velocity u and 
for the magnetic field B may be written in the following 
form, valid at time scales 3> fJ" 1 (f2j = eB/rriiC is the ion 
cyclotron frequency) and spatial scales 3> Pi — Vth,i/^i, 



du 

dB 

~dt 





B 2 \ 




fP-L 




+ V- 



bb (p ± - p|| + B 2 



B Vu BV u, 



,(1) 
(2) 



where d/dt = d/dt + u-'V is the convective derivative, p± 
and p» are plasma pressures perpendicular and parallel to 
the local direction of the magnetic field B, respectively, 
b = B/B, the factor of has been absorbed into 

B, and the resistive term has been omitted in Eq. J3J in 
view of the tiny value of the resistivity. The turbulent 
motions in clusters are subsonic (U < i>th,i)j so we may 
take V • u = and set p = 1. The magnetic field is, thus, 
in units of velocity, pressure in units of velocity squared. 

The proper way to compute p± and p» is by a kinetic 
calculation. In the collisional limit, this was done in Bra- 
ginskii's classic paperi22i It is instructive to obtain his 
result in the following heuristic way that highlights the 
physics behind the formalism^ Charged particles mov- 
ing in a magnetic field conserve their first adiabatic in- 
variant p = niiv\/2B. When A m f p S> Pi, this conserva- 
tion is only weakly broken by collisions. As long as p is 
conserved, any change in the field strength causes a pro- 
portional change in p±: summing up the first adiabatic 
invariants of all particles, we get px/B = const. Then 



1 dpj 



P± 



dt 



1 dB 
~B~dl 



— —r. Vi, 



Pi- ~P\\ 
P± 



(3) 



where the second term on the right-hand sight represents 
the relaxation of the pressure anisotropy at the ion colli- 
sion rate va ~ Uth.j/Amfp^" Using Eq. @ with V • u = 
and balancing the terms in the rhs of Eq. (0, we get 



P± ~P\\ ~ M| 



1 dB 
B~dT 



p iibb : Vu, 



(4) 



where fiu ~ px/^u ~ i'th,iA m f p is the "parallel viscos- 
ity." Equation turns out to be quantitatively correct 
provided the value of /zy with the correct prefactor cal- 
culated by Braginskii22 is used. Thus, the emergence of 
the pressure anisotropy is a natural consequence of the 
changes in the magnetic-field strength and vice versaS 
The energy conservation law based on Eqs. with 
V • u = and on Eq. Q is 



d_ 

dt 



(u 2 ) , (B 2 ) 



= -M||<|bb: Vu| 2 ) 




(5) 



Thus, the Braginskii viscosity only dissipates such mo- 
tions that change the strength of the magnetic field. Mo- 
tions that do not affect B are allowed at subviscous scales. 
In the weak-field regime, these motions take the form of 
plasma instabilities. When the magnetic field is strong, a 
cascade of Alfven waves can be set up below the viscous 
scale. Let us elaborate. 

The simplest way to see that pressure anisotropies lead 
to instabilities is as follows Imagine that the large-scale 
energy sources stir up a "fluid" turbulence with u, p±, 
pn, B at time and spatial scales above viscous. Would 
such a solution be stable with respect to much higher- 
frequency and smaller-scale perturbations? Linearizing 
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Eq. and denoting perturbations by 5, we get 
- iujSu = -ik (8 P x + BSB) + (p x - p l{ + B 2 ) SK 



+ ibk 



Sp ± - 5p\\ - (px - P\\ - B 2 ) 



5B 
~B 



,(6) 



where K = b • Vb is the curvature of the field. Using 
Eq. J2J, SK = k 2 Su±/iuj. We can immediately separate 

the Alfven- wave-polarized perturbations (Su cx k x b), 
for which we obtain the dispersion relation 



±/c|| (p x -P\\+B 2 



1/2 



(I) 



When p|| — p± > B 2 , u> is purely imaginary and we have 
what is known as the firehose instabilityiSiSiSiSi The 
growth rate of the instability is oc k\\ , which means that 
the fastest-growing perturbations are at scales far be- 
low the viscous scale or the mean free path. Therefore, 
Eqs. (^[3} are iU posed wherever pu — p± > B 2 . To 
take into account the instability and its impact on the 
large-scale dynamics, a kinetic, rather than fluid, descrip- 
tion must be adopted. A linear calculation based on the 
hot plasma dispersion relation shows that the instabil- 
ity growth rate peaks at the ion gyroscale, fey p.; ~ 1. 
While the firehose instability occurs in regions where the 
magnetic- field strength is decreasing [Eq. (0}], a kinetic 
calculation of Spx and Sp\\ in Eq. © yields another insta- 
bility, called the mirror mode^ that is triggered wher- 
ever the field is increasing (px > pu). Its growth rate 
also peaks at kup% ~ 1. 

To sum up briefly, external energy sources drive ran- 
dom motions, which change the field [Eq. J2J], which 
gives rise to pressure anisotropics [Eq. 10}], which trig- 
ger the instabilities. The latter are stabilized when 
B 2 > \p± — pu | and the firehose fluctuations, in particu- 
lar, mutate into Alfven waves, which can cascade with- 
out being affected by collisions all the way to the ion 
gyroscale We shall revisit the subject of the Alfvcn- 
wave cascade in Sec. IVII but first let us see what the 
extraordinarily unstable nature of the ICM in the weak- 
field regime implies. 

What it implies is, in fact, not altogether clear, as the 
quantitative theory of the evolution of the instabilities 
beyond the linear stage is a difficult taskiSi While the 
standard quasilincar scheme can be implemented more 
or less rigorously in the limit of small anisotropy for par- 
ticular cases in which the instability growth rates peak 
at scales much larger than the ion gyroscale (specifically, 
when k± = for the firehose and fcj_ 3> fcii for the 
mirror)^ it is not clear that this is relevant, because in 
the general case, the growth rates peak at kpi ~ 1. The 
difference between this situation and the case of kpi <C 1 
is qualitative: for the fluctuations at the gyroscale, the p 
conservation is broken, so they can give rise to effective 
particle scattering, while the quasilincar saturation of the 
fluctuations with kupi <C 1 happens essentially because 
Spx and 5p\\ offset the initial pressure anisotropy^ 




-10 -8 -6 -4 

lg(B/B eq ) 

FIG. 1: Re c ff(B): solution of Eq. (f 1 1 pi calculated numerically 
for the cool-core parameters of Table Q] 



Here we would like to leapfrog the rather forbidding 
task of constructing a quasilincar theory based on the 
hot plasma dispersion relation with kpi ~ 1 and in- 
stead simply assume that the istabilitics saturate with 
(SB 2 )/ B 2 ~ (|A| - 2/(3) ai , where A = (px ~ P\\)/px, 
(3 = 2px/B 2 , and ot\ is some positive power. We further 
assume that the effect of the saturated instabilities is 
to scatter particles at the rate Scatter ~ lmax(SB 2 )/B 2 , 
where 7 max ~ (|A| — 2//3J 2 is the maximum growth 
rate of the instabilities (a.2 = 1 for the mirror and 
a- 2 = 1/2 for the firehose [Eq. @]). We may then con- 
struct the effective collision rate 



^eff ~ V%i + ^scatter 



A| 



(8) 



where a — ol\ + ct2 > 0. This changes the characteristics 
of the turbulence: the effective mean free path of the 
particles is A m f PiG ff ~ vth.i/veS, the effective (parallel) 
viscosity of the ICM is p\\ te s ~ ^th.i^mfp.eff and, therefore, 
the effective Reynolds number is 



Re, 



eft 



UL 

M||,eff 



UL 

—o feff- 



(9) 



On the other hand, using Eq. with effective viscosity 



1/2 

/i|| jeff and |Vu| - (U/L)Re c ' s , we get 

( U \* 1/2 



Rel 



(10) 



From Eqs. l(8T41()|l. we assemble a (model) equation for 
the effective Reynolds number: 



Re 



eff 



Re 



L 



Pi,eq \Vth,i 



U 



2a + l 



1 



VRccff 



2B 2 . 

) — .(n) 



B 
B, 



CC] 



where Re is the original Reynolds number based on col- 
lisions and Pi^cq is the ion gyroradius for B = B cq . The 
assumption that Re c ff adjusts to the value of B instanta- 
neously is justified by the fact that the instability growth 
rates are much faster than all other relevant time scales. 
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FIG. 2: The mechanism of the fluctuation dynamo. 



Equation models qualitatively the assumed out- 
come of an as yet inexistent proper theory of the viscos- 
ity of magnetized ICM. Its solution is plotted in Fig. ^ 
We shall use this model shortly in our discussion of the 
fluctuation dynamo in clusters. 



V. FLUCTUATION DYNAMO 

In Sec. IIIII we reviewed the observational evidence that 
testified to the presence of a dynamically significant ran- 
domly tangled magnetic field in clusters. What is the 
origin of this field? There arc numerous physical rea- 
sons to expect that a certain amount of seed magnetic 
energy predates structure formation and was, therefore, 
already present at the birth of clusters i 4 ^*^ Typical val- 
ues given for the strength of such field are in the range 
of B sccd ~ 1(T 21 - lCT 17 , although this may be an 
underestimated It then falls to the random motions of 
the cluster plasma to amplify the field to its observed 
magnitude of a few /iGi££ This, indeed, they should be 
able to do by means of the fluctuation (or small-scale) 
dynamo mechanism: the random stretching of the field. 
It is a fundamental property of a succession of random (in 
time) linear shears that it leads on the average to expo- 
nential growth of the energy of the magnetic field frozen 
into the mcdium4244i2£ The rate of growth is roughly 
equal to the rate of strain (shear, or stretching rate) of 
the random flow. While the mathematical theory of this 
process can be nontrivial^ the physics of it is basically 
illustrated by Fig. 

In Kolmogorov turbulence, the rate of strain is domi- 
nated by the viscous scale, so | Vu| ~ t^ c ~ ([//L)Re 1 ^ 2 . 
In fact, what is relevant for the growth of the magnetic 
field is not the full rate-of-strain tensor but its "paral- 
lel" component, bb : Vu. Since this is exactly the type 
of motion damped by Braginskii viscosity [see Eqs. 
EJ], we can, for the purposes of the fluctuation dynamo, 
ignore any subviscous-scale velocity fluctuations. Thus, 
the magnetic field should grow according toS 



1 dB 
~B~dt 



bb : Vu 



U 

T 



Re 1 / 2 . 



[12) 



In order to be useful, this amplification has to be done 
over a time that is shorter than the typical age of the 
clusters, i.e., a few Gyr. As L/U ~ 1 Gyr and as 
the Reynolds number based on collisional ICM viscos- 
ity is quite low, it has been a concern of many au- 



thors that the cluster life time might not be (or, pos- 
sibly, is only just) sufficient to bring the field to its 
current s trengthi2iiIi22ii2i52i5LS (this obviously depends 
on the magnitude of the seed field). Similar issues 
arise in simulations of the growth of magnetic fields in 
clusters! ^ 4 ^'^? However, if Eq. (fTTJ) is a good model 
of the ICM viscosity, the value of the Reynolds number 
based on the collisional viscosity has to be amended when 
the ICM is magnetized: since the prefactor L/pi teq in the 
second, magnetic-field dependent term in Eq. Ijllj) is ex- 
tremely large (see Table HJ , this term will dominate al- 
ready for very small values of the magnetic field strength, 
namely, for B > B\, where 



Bi = B, 



cq 



Pi,eq ( ^th,i 



l+2a 



Re 



l+a/2 



(13) 



B\ is typically so small (Table HJ) that the seed field 
in clusters may already exceed this value. Thus, the 
effective Reynolds number of the ICM will be much 
larger than the one based on collisions. Substituting 
Re c ff instead of Re into Eq. (|12fl . we see that, since 
Re e ff oc i3 2 /( 2 +") in the weak-field regime, the dynamo is 
sclf-accclcrating and the field will increase not exponen- 
tially but explosively: 



B(t) 



B(0) 



(1-t/tc) 



2+Q ' 



(14) 



where B(0) ~ the greater of -B S eed and B\ and t c — 
(2 + a)(L/[/)Re- 1/2 [5i/5(0)] 1/(2+Q) is at most (for 
Bseed < B\) the viscous turnover time i v isc associated 
with the collisional Re. The explosive stage contin- 
ues until the amplified field starts suppressing the in- 
stabilities, i.e., when (3 drops to values comparable to 
{vth,i/U) 2 Kel(^ . This happens at B ~ £?2, where 



B, = R 



cq 



. l+2a 
Vtb,i \ Pi,eq 



l/(5+2a) 



(15) 



Thus, we have a mechanism that amplifies the field by 
many orders of magnitude from any strength above B\ to 
Bi in finite, cosmologically short, time ~ t c . The value of 
Bi turns out to be only just over an order of magnitude 
below B eq (see Table QJ. 

Further growth of the field is algebraically slow (B ~ 
i 1 / 2 ), but it does not have to go on for a very long 
time because B^ is already quite close to the observed 
field strength. To be precise, there are two algebraic 
regimes. During the first, Re e ff is still controlled by the 
second term in Eq. (|llf) as B hovers just below £? cq Re^/ 4 
while B is increasing and Re c ff is decreasing. Eventually, 
B ~ Bvisc = -B cq Re -1 / 4 and Re ff is returned back to Re 
(plasma instabilities are suppressed) J£ As this is also the 
field strength at which the field has energy comparable 
to the energy of the viscous-scale motions, any further 
growth of the field is a nonlinear process, in which the 
back reaction of the field on the flow has to be taken into 
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FIG. 3: Evolution of the magnetic-field strength for the cool- 
core parameters of Table [I] 




FIG. 4: Cross sections of |u| (left) and |B| (right) in the 
saturated state of a simulation with Re ~ 100, Rm = 1000 
(run B of Ref.|H). 



account. This can be done by assuming that, as the field 
grows, it can no longer be stretched by motions whose 
energy it exceeds and that, therefore, at any given time, 
the dominant stretching is exerted by motions whose en- 
ergy is equal to that of the field. Denoting their velocity 
and scale by ui and I and using uf ~ B 2 , we haveSS 



dt 



-B 



11 
I 



B 2 



I 



const, 



(16) 



where e is the energy flux through scale I, which, by 
Kolmogorov's assumption, is independent of I. Thus, 
B ~ (et) 1 / 2 , the second algebraic regime. Obviously, 
this can only go on until I ~ L, B ~ B cq , and no further 
amplification is possible. 

The time evolution of the field is illustrated in Fig. [3J 
where we plot the result of a numerical integration of 
Eq. (|12|) with Re replaced by Re ff. To implement the 
idea of the nonlinear suppression of stretching by motions 
whose energy is smaller than B 2 , we have amended the 
definition of Re e ff in the following simple wayS 



Re 



(nonlin) 
eff 



Re of t 



I- A 



1 



[1 + Re 



1/2 
efit 



B 2 /B 2 q ) 2 



A 



(17) 



where A = 1/(1 
when B — > and Re 



Re^ 2 ) 2 , so that Reiff° nl = Re ff 



eft ; ' DU ullaL ilc eff 
fnonlm) _^ q w \ lcn Q _> £ gq 

determined for each value of B by Eq. pijl. 



Re e ff is 



VI. MAGNETIZED PLASMA TURBULENCE 

Thus, the clusters should have no trouble developing 
magnetic fields of observed strength over just a fraction 
of their life timeS This, however, is only one part of the 
problem. The other, much more difficult, part is to deter- 
mine the spatial structure of the field. ThcorySfliSi and 
numerical simulations^iS of the fluctuation dynamo in 
a conventional MHD formulation with fixed isotropic vis- 
cosity and much smaller resistivity (Rm 3> Re) produce 
magnetic fields arranged in folded flux sheets whose par- 
allel (locally to themselves) scale is similar to the scale of 



the velocity that does the stretching, while the field re- 
versal scale in the direction perpendicular to itself is the 
resistive scale. Fig. illustrates why such a folded struc- 
ture must be a generic outcome of random stretching. 
Less obviously, it turns outS that the folded structure is 
a feature not only of the kinematic growth stage of the 
dynamo but also of its saturated state (Fig.^J. In the lat- 
ter case, the fold length is, thus, l» ~ L, while the rever- 

1/2 

sal scale is lj_ ~ l rcs ~ LRm ' . The magnetic-energy 
spectra of such fields appear to peak around k ~ 1/Zj_. 

This prediction is clearly not in agreement with what 
seems to be the observationally supported picture of the 
magnetic fields in clusters. While the typical scale Ib of 
the cluster fields is significantly smaller than the outer 
scale L of the turbulent velocities, it is certainly not the 
resistive scale — at least not the one based on the stan- 
dard Spitzcr resistivity. What then determines the rever- 
sal scale in clusters? While the definite solution of this 
puzzle remains elusive, we offer the following qualitative 
argument that at least constrains the answer. 

Examining again Fig. [21 it should be obvious that if 
the folded magnetic field is to reverse its direction, it 
must turn a corner somewhere and that in that corner, 
the field should be weaker than in the straight segments 
of the fold. Furthermore, as there arc antiparallcl fields 
in the straight segments themselves, there must be layers 
of weak field between them (assuming that folds are flux 
sheets: i.e., that the antiparallel fields in Fig. [21 are, in- 
deed, approximately aligned). These are, in fact, also of 
the "corner" type: this becomes clear if the fold depicted 
in Fig. [21 is thought of as a fragment of a larger nested 
folded structure formed by repeated stretching/shearing 
and bending. This view of the dynamo-generated fields 
as either strong and straight or weak and curved is con- 
firmed in numerical simulations, which find an almost 
perfect anticorrelation between the field's strength and 
its curvature^ While the exact quantitative relation be- 
tween the field strength and its scale may be a non- 
trivial issue, here we shall continue in the spirit of in- 
tuitive reasoning and argue that by flux conservation, 

£ s traight/-B corn cr ~ hl^ 1 - ( tnis a S am re hes on assuming 
that folds are flux sheets). Thus, the larger is the aspect 
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FIG. 5: A schematic illustration of the structure of cluster 
turbulence proposed in Sec. I VII 

ratio l\\/l±, the larger is the contrast between the strong 
field in the straight segments of the folds and the weak 
field in the corners. We assume that the rms value of 
the field is determined by the straight segments because 
these are the fields that are stretched by turbulence. It is 
their growth that we studied in Sec. [3 In the saturated 
state, we expect -B s traight ~ B eq . For such a strong field, 
the plasma instabilities are suppressed. It is intuitively 
clear that they must be suppressed in the regions of the 
weak field as well, i.e., the field there cannot be weaker 
than i?2 [Eq. (jl5|l ]. Indeed, as we saw in Sec. E| the 
explosive dynamo mechanism brings any field up to this 
value nearly instantaneously (for B > i?2, the growth 
is much slower). We may conjecture that the maximum 
aspect ratio of the folds is set by the maximum contrast 
in the field strength l\\/l±_ ~ B eq /B 2 ~ Substituting the 
numbers, we find (Table P) that this prediction gives the 
field reversal scale no smaller than a few per cent of the 
outer scale L (taking lu ~ L). This is in passable agree- 
ment with the observational evidence, which is the best 
one can expect, given the highly imprecise nature of our 
argument, of most observational inferences, and of the 
definitions of such quantities as Ib, h and L. For 
comparison of our model with observational data for a 
number of individual clusters, see Ref. l20i 

It is fair to acknowledge that the above argument, 
while providing a useful constraint, falls short of a sat- 
isfactory explanation of the field structure. One might 
argue that if, in the course of the turbulent stretch- 
ing/shearing of the field, a region of field strength below 
Bi appears (as explained above, in the corner of a fold), 
Re e ff there becomes very large and a localized spot of 
high-Reynolds-number turbulence is formed. This should 
have two principal effects. The first is akin to that of a 
locally enhanced turbulent resistivity, so the field that 
violates our constraint is continuously destroyed. The 
second is a burst of explosive fluctuation dynamo in the 
spot, which produces more folded field with B > Bi 
and thus shuts itself down. These folds arc then further 
stretched, sheared, etc., again subject to the constraint 
that they are destroyed and replaced by new ones wher- 
ever a spot of weak field appears. We do not currently 
have a more detailed mechanistic scenario of how exactly 



the folded structure with l\\/l± ~ B eq /B2 is established. 
It may be feasible to test these ideas numerically by solv- 
ing MHD equations with viscosity locally determined by 
the magnetic-field strength according to Eq. (|11|) . 

Let us assume that cluster fields do indeed have a 
folded structure with a direction reversal scale l± ~ 
(0.01...0.1)L, possibly determined by the argument given 
above^ The magnetic-energy spectrum then peaks at 
k ~ What is the structure of the turbulence above 

and below this scale? At scales I <C l±, the magnetic 
field reversing at the scale l± will appear uniform and, in 
accordance with the old idea of KraichnanS could sup- 
port a cascade of Alfven waves. This cascade can rig- 
orously be shown to be described by the equations of 
Reduced MHD at collisionless scales all the way down 
to the ion gyroscale^ The currently accepted theory of 
such a cascade, primarily associated with the names of 
Goldreich and Sridhar^ is based on the conjecture that 
at each scale, the Alfven frequency is equal to the turbu- 
lent decorrelation rate. The result is a fc -5 / 3 spectrum of 
Alfvenic fluctuations — this possibly explains what ap- 
pears to be a A; -5 / 3 tail in the observed spectrum of mag- 
netic energy for the Hydra A corei We cannot embark 
on a detailed discussion of the theory of the Alfven-wave 
cascade here, so the reader is referred to Ref. |6^ for a re- 
view and to Ref. |3l] for the theoretical basis of extending 
this theory to collisionless scales. 

Above the reversal scale, I 3> l±, the cluster turbu- 
lence should resemble the saturated state of isotropic 
MHD turbulence: a magnetic-energy spectrum with a 
positive spectral index corresponding to folded fields and 
a kinetic-energy spectrum populated in the incrtial range 
by a peculiar type of Alfven waves that propagate along 
the folds (i.e., simultaneously perturbing the antiparallel 
magnetic field lines) This type of turbulence is also 
reviewed in Ref. 16a . It is probably of limited relevance 
for clusters because the Reynolds number in the ICM is 
not large enough to allow a well-developed incrtial range. 

Fig. [S] summarizes the — admittedly, rather specula- 
tive — picture of cluster turbulence proposed above. We 
offer this sketch in lieu of conclusions. While we believe 
that the set of physical arguments that has led to it is not 
without merit, it is clear that much analytical, numeri- 
cal and observational work is needed before a conclusion 
can truly be reached in the study of turbulence, magnetic 
fields and plasma physics in clusters of galaxies. 
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The cool cores occupy roughly the central 100 kpc, com- 
pared with the overall cluster size of order 1 Mpc. Such 
a structure is, in fact, characteristic of the older type of 
clusters known historically as cooling-flow clusters (the 
more observationally popular members of this group are 
Hydra A, Centaurus and Perseus), while younger, less dy- 
namically relaxed clusters (e.g., Coma) have flatter density 
and temperature profies. A good summary of the parame- 
ters in a number of cooling-flow clusters and the relevant 
references can be found in Ref. l20t 

These doubts are reinforced by the similarity between the 
inferred flow patterns associated with rising bubbles and 
known such patterns in viscous, rather than turbulent, 
fluid flows-ifi 

An estimate of resistivity n based on the standard Spitzer 
formula leads to the enormous values of the magnetic 

Reynolds number, Rm ~ UL/n, g iven in Table Q 

http://www.lofar.org/ ; http://www.skatelescope.org/ 
This is only valid if the characteristic parallel scales fcjT of 
all fields are larger than A m f p . In the collisionless regime, 
k\\\m{ P S> 1, we may assume that the pressure anisotropy 
is relaxed in the time particles streaming along the field 
cover the distance k^ 1 : this entails replacing vu in Eq. © 
by fc||«th,i- 

The anisotropy is small: A = (p± — p\\)/p± ~ 
(U/v t h,i) 2 Re~ 1/2 . It turns out that (f7/«th,i)Re~ 1/4 is the 
natural small parameter that can be used to develop a re- 
duced kinetic theory for the cluster plasmap— 
In space physics, where the importance of the plasma in- 
stabilities driven by pressure anisotropies of what in their 
case is an almost entirely collisionless plasma has long been 
understoodi— A vast but inconclusive literature exists on 
this subject. 



Several authors have argued that no amplification is, in 
fact, necessary either because the seed fields may have al- 
ready been strong enough to account, after compression 
of the cosmological plasma into cluster, for the observed 
magnitude of the cluster field^- or because the field could 
be generated in AGNs and then ejected into the ICMp— 
We forego the discussion of these possibilities, primarily 
because we find the idea that the ICM turbulence will pro- 
duce the right amount of magnetic energy independently 
of its history or external circumstances more appealing on 
a fundamental physics level and, indeed, intuitively sup- 
ported by the fact that the energy density of the observed 
field is close to the energy density of the fluid motions. 
Note that turbulence in the sense of a broad inertial range 
is not required for the fluctuation dynamo. The viscous- 
scale motions, which dominantly stretch the field, are 
random but spatially smooth. Whether the viscous scale 
ivisc ~ I/Re _3//4 is much smaller than the outer scale L 
is inessetial as long as the motions are random — not a 
problem in clusters, given the vigorous random stirring. 
Thus, the Reynolds number in all of our calculations need 
not be large. Numerical simulations of a randomly stirred 
magnetohydrodynamic (MHD) fluid with Re in the range 
1 — 10' 5 confirm the insensitivity of the dynamo effect to 
the value of Ra^S. 

This, effectively, is the regime assumed by Sharma et aim— 
in their numerical simulations of the magnetorotational 
instability in a collisionless plasma. Their set of fluidlike 
equations incorporates anisotropic pressure and an effec- 
tive collision rate that enforces the marginal state of the 
plasma instabilities: (v t h,i/U) 2 Re^ ~ (3 in our notation. 
This means, in particular, that young clusters should al- 
ready have dynamically significant magnetic fields. 
This does not mean that cluster fields must be volume fill- 
ing. Indeed, a real cluster is probably a patchwork of tur- 
bulent and quiescent regions (or rather regions with widely 
varying rms rates of strain) rather than a volume homoge- 
neously filled with turbulence. Furthermore, the turbulent 
dynamo only produces magnetic fields everywhere on the 
average (in time), while any particular snapshot is very 
intermittenli— — see, e.g., Fig. 2] 

Just such a structure appears to be evidenced by the po- 
larised emission map from a radio relic in cluster A2256p— 



